Forward scattering amplitudes and the thermal operator representation 
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We develop systematically to all orders the forward scattering description for retarded amplitudes 
in field theories at zero temperature. Subsequently, through the application of the thermal opera- 
tor, we establish the forward scattering description at finite temperature. We argue that, beyond 
providing a graphical relation between the zero temperature and the finite temperature amplitudes, 
this method is calculationally quite useful. As an example, we derive the important features of the 
one loop retarded gluon self-energy in the hard thermal loop approximation from the corresponding 
properties of the zero temperature amplitude. 
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I. INTRODUCTION 

In a series of papers the idea of a thermal 

operator representation 0, |g has been developed exten- 
sively in both the imaginary time formalism as well as the 
real time formalism of closed time path. In simple terms, 
the thermal operator representation relates a Feynman 
graph at finite temperature (with or without a chemical 
potential) to the corresponding graph at zero temper- 
ature. As we have argued earlier, the thermal operator 
representation offers a powerful method for studying var- 
ious questions at finite temperature. As an example, we 
have shown in an earlier paper "t^I how the cutting rules 
at finite temperature (with or without a chemical poten- 
tial), in the closed time path formalism, can be derived 
starting from those at zero temperature. This derivation 
also clarifies the miraculous cancellations that arise in 
an explicit demonstration of a cutting description for the 
imaginary part of a thermal amplitude |3, |9| • 

At finite temperature, retarded amplitudes play a sig- 
nificant role in studying various physical phenomena. 
Plasma oscillations provide a very simple example of this. 
When a thermal plasma is perturbed weakly, the sub- 
sequent response of the plasma to the pert urbation is 
studied using the linear response theory \WL UlL H^ . In 
particular, the damping of the oscillation in the plasma is 
understood by analyzing the poles of the retarded prop- 
agators of the particles moving through the plasma. Of 
course, at finite temperature very few quantities can be 
evaluated exactly, but the forms of thermal amplitudes 
simplify considerably either in the low temperature or the 
high temperature limits. In many phenomena of physical 
interest (such as quark-gluon plasma phase transitions, 
early universe etc), it is the high temperature behavior 
that is relevant. While there are many ways of evaluating 
the high temperature behavior (also known as the hard 
thermal loop approximation 13] ) of thermal amplitudes, 
the forward scattering description for the retarded ampli- 
tudes provide an efficient calculational tool 0| . This can 
be seen from the following simple example. Let us con- 
sider a scalar field theory with a cubic interaction in six 



dimensions (which is similar to non-Abelian gauge theo- 
ries in four dimensions). The thermal correction to the 
one loop retarded self-energy can be directly calculated 
(see Fig. 1) in this theory and, after doing the internal 
energy integral, leads to 
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FIG. 1: One-loop retarded self-energy in 0^ theory. 
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thermore, nB{Ek) denotes the bosonic distribution func- 
tion and po is assumed to correspond to po -\- ie which 
is necessary for the retarded self-energy. At very high 
temperatures where \k\ ^ p^,m, the masses can be ne- 
glected and we see from Q that the high temperature 
limit needs to be calculated carefully since there are en- 
ergy differences in the denominator. 

On the other hand, the forward scattering description 
of the same retarded self-energy involves diagrams where 
one of the internal propagators in the loop is thermal and 
on-shell while the other corresponds to a zero tempera- 
ture retarded propagator. They are shown in Fig. 2 and 
lead immediately to 
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FIG. 2: Two forward scattering amplitudes for the one-loop 
retarded self-energy. 
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Here the "ze" in po (denoting a retarded propagator) as 
well as fco = are to be understood. At high tempera- 
ture where masses can be neglected, this takes the simple 
form 
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The structure resulting directly from the forward scatter- 
ing description is very simple and interesting (of course, 
the same structure would also arise in a direct calcula- 
tion, but limits have to be taken carefully and terms have 
to be grouped properly before this simple structure is ob- 
tained). First of all, we note that the coefficient of the 

term "°^^J^^'* in the integrand is manifestly Lorentz covari- 

ant and is a homogeneous function of degree zero in the 
external momentum and of degree (—2) in the internal 
momentum. The manifest Lorentz covariance is broken 
at high temperature only when the angular integration is 
carried out. In fact, if we carry out the integration over 
|fc|, the high temperature limit is obtained to be 
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where we have defined fc^ = (l,fc). This Lorentz co- 
variant structure of the integrand in (which results 
directly in the forward scattering description) is very 
helpful and has been used to derive in a simple way, in 
the hard thermal loop approximation, the effective ac- 
tion for QCD as well as the energy-momentum tensor for 
the quark-gluon plasma . This method is also conve- 
nient for the analysis of the high temperature behavior 
of gauge field theories in a curved space-time . It is 
also worth noting that the study of the solution of the 
transport equation at high temperature leads to struc- 
tures naturally arising in the forward scattering method 

In spite of its success, a simple and general derivation 
of the forward scattering amplitudes to all orders at finite 



temperature is so far lacking. As we have already argued, 
the thermal operator representation P, 0, provides a 
powerful method for obtaining results at finite tempera- 
ture starting from zero temperature and in this paper we 
show how the forward scattering amplitudes for retarded 
thermal n-point functions can be derived through the use 
of the thermal operator representation. The thermal op- 
erator representation is clearly meaningful in studying 
this question if there exists a forward scattering descrip- 
tion at zero temperature. In this paper we derive the 
forward scattering description for retarded amplitudes in 
zero temperature field theories. The thermal operator 
representation then directly leads to the forward scat- 
tering description at finite temperature and clarifies the 
origin of the nice structures observed in the context of 
the forward scattering amplitudes at high temperature. 

The paper is organized as follows. In section II, we de- 
velop the forward scattering description for retarded am- 
plitudes of a scalar field theory at zero temperature. In 
section III, we show how the thermal operator represen- 
tation leads directly to the forward scattering description 
for retarded thermal amplitudes. In this section, we also 
point out various interesting features of retarded ampli- 
tudes both at zero as well as at finite temperature. In 
section IV we discuss the forward scattering description 
for the Yang-Mills theory. In particular, we emphasize 
that various nice properties in the hard thermal loop ap- 
proximation such as transversality, manifest Lorentz co- 
variance and gauge invariance of the integrand of the one 
loop retarded self-energy follow simply from the proper- 
ties of the zero temperature amplitude through the ther- 
mal operator representation. We conclude with a brief 
summary in section V. 



II. FORWARD SCATTERING DESCRIPTION 
AT ZERO TEMPERATURE 

The idea of forward scattering description basically al- 
ready exists even at zero temperature although it is not as 
well developed and is certainly not widely known. There- 
fore, in this section, we will develop the idea of forward 
scattering amplitudes at zero temperature systematically 
for retarded amplitudes so that the thermal operator rep- 
resentation can lead directly to the forward scattering 
amplitudes at finite temperature. The basic idea behind 
a forward scattering description at zero temperature 
is the simple observation that a (time ordered) Feynman 
propagator (for a massive scalar particle, for simplicity) 
can be expressed as 

2Tr0{~ko)S{k^ -m^), 



+ ie (fco + ie)2 
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where we have identified E^. = k^ 



(5) 

m^. Namely, the 
Feynman propagator is the sum of the retarded propa- 
gator and a negative energy propagator. As a result, if 
we have a simple one loop diagram with n scalar propa- 
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gators, the amplitude (neglecting vertex factors) can be written as 
I 
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where we have denoted the momentum in the zth propa- 
gator as ki which is a sum of the loop momentum k and 
some combination of the external momenta whose ex- 
plicit form is not relevant for our discussion. The form of 
the integrand in ((BJ is quite interesting. The first term 
in the product which will involve only products of re- 
tarded propagators would vanish when integrated over 
energy (which can be seen simply as a consequence of 
the fact that all the poles lie on the lower half of the 
complex plane and, therefore, the contour can be closed 
in the upper half plane to yield zero). The other terms in 
the expansion of the right hand side would involve terms 
with a number of retarded propagators and the rest of 
the propagators on-shcU. If we assume that an on-shell 
propagator can be thought of as a cut open line repre- 
senting an on-shell particle coming in and going out, this 
is very roughly a forward scattering description, namely, 
a Feynman amplitude can be expressed as a sum of di- 
agrams that involve a number of on-shell particles scat- 
tering in the forward direction (their momenta are un- 
changed in the scattering) and retarded propagators. It 
is worth noting from the form of JHJl that the series of 
forward scattering diagrams may involve completely dis- 
connected diagrams (which is not the case for retarded 
amplitudes at finite temperature), but we would like to 
point out that this is only a consequence of the fact that 
we are looking at a time ordered Feynman amplitude. 

On the other hand, we are interested in retarded am- 
plitudes and, as is well known, these are quite hard to 
construct at zero temperature within the context of the 
conventional Feynman propagators. However, if we dou- 
ble the degrees of freedom (the theory with the doubled 
degrees of freedom can be taken as the zero temperature 
limit of the theory at finite temperature in the closed 
time path formalism as described in [3)@])j a- diagram- 
matic representation of retarded amplitudes can be con- 
structed in a straight forward manner. With this in mind, 
let us look at a scalar field theory with a 0'^ interaction 
with doubled degrees of freedom and we denote the two 
field degrees of freedom as 0+ and 0_ . The propagator 
for the doubled theory corresponds to a 2 x 2 matrix 



A = 



(7) 



and in the momentum space, the components take the 
forms 
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A+_(p) = 2TTe{~p^)5[p' -m'), 
A_+(p) = 2^0(po)<5(p' - ™2), 

A (p) = lim — 



(8) 



(Unlike in our earlier papers 0, |3i j here we will follow 
the simple convention of representing quantities at zero 
temperature without any superscript (T = 0). We will 
reserve the superscript (T) only for quantities at nonzero 
temperature in order to simplify the notation.) 

However, since the forward scattering amplitudes have 
a physical description in the mixed space (we will discuss 
this later), we will analyze the problem in this context 
(the momentum space analysis that is normally done can 
be obtained from our results through a Fourier trans- 
formation). In the mixed space, the components of the 
propagator can be obtained from a Fourier transforma- 
tion of © and have the forms (see also for various 
notations and conventions) 
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where E — \/p^ + "rn^ describes the on-shell energy of 
the particle and we are suppressing the "ie" in the expo- 
nents for simplicity. Thus, we see that A±zp describe re- 
spectively the on-shell negative and positive energy prop- 
agators. The vertices involving the <j)^ fields and the 
fields differ by a relative negative sign. 

The time ordered components of the propagators in @ 
satisfy the constraint 



A__ = A, _ + A_ 
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(10) 
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It is now simple to see that the retarded and the advanced 
propagators of the theory can be identified with 

- 0{t) ^ (e-^* - e^^*) , 
A!,{t,E) = A++{t,E)-A^+{t,E) 

= 0{-t)^{e''''-e-^'='), (11) 

much hke at finite temperature (Qj. From these defini- 
tions, we note that 

An{~t,E) = AA{t,E), (12) 

as we would expect. As a result of these relations, we 
can decompose the matrix propagator in Q as 

A = P + A+_Q, (13) 

where 

= (;;). (14) 

I 
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FIG. 3: Sum of two diagrams which gr 



Let us next note that a retarded n-point amplitude at 
any loop can be defined as follows. If we assume that the 
time ti corresponding to the first index of the amplitude 
is the largest among the time coordinates, then we have 

^n,R{tl, ■ ■ ■ ,tn) — ^ F+ai---a„_i(il, • ■ • , ^n) (15) 

o,=± 

where we have suppressed the energy dependence of the 
amplitude for simplicity. Here at denote the "thermal 
indices" of the fields which can take the values "±" . For 
example, the retarded two point function (self-energy) at 
one loop would correspond to the sum of the two dia- 
grams in Fig. 3. Let us also note here for future use 
that, for any n-point amplitude, 

^ ^ r(jj^(j203 - a„ — Oj (16) 
o. = ± 

which follows from the largest time equation We are 
now ready to derive the forward scattering description 
for retarded amplitudes to all orders at zero temperature 
and we do so in two steps. 



k. El 




k + P,E2 



the retarded self-energy at one-loop. 



A. Forward scattering amplitudes at one loop 

The forward scattering amplitudes for retarded n-point 
functions can be derived algebraically at one loop (which 
is the reason for separating the derivation into two cases). 
First we note from H15(l that the retarded amplitude con- 
sists of terms where each vertex other than the largest 
time is summed over the thermal index "±" . Further- 
more, as we have already pointed out, the vertex for the 
(j>- field has a relative negative sign compared to that for 
the 0+ field. Thus, summing over the thermal index of 
the graph at one loop can be effected by multiplying the 
matrix propagator with a 2 x 2 matrix at the vertex 
where the thermal index is being summed. For exam- 
ple, for the case of the retarded self energy at one loop 



(see Fig. 3), we note that (once again we are neglecting 
factors associated with the vertices as well as the depen- 
dence on external energies for simplicity) 

4'^(ii-t2) 

[A++{ti - t2, Sl)A++(t2 - ti,E2) 

-A+^{h-t2,E,)A^+{t2-ti,E2)] 

^[A(ii -t2,i?l)<T3A(t2 -tl,i?2)]++.(17) 

Note also that the retarded amplitude is obtained by tak- 
ing the "H — h" component in the matrix product (simply 
because we start from a vertex and end at the same 
vertex). As a result of this simplification, the retarded 
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n-point amplitude at one loop, shown in Fig. 4, can be 
written as 




FIG. 4: One-loop diagram for the retarded n-point function. 
Sum over the thermal indices from t2 to t„ is to be understood. 
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We can use the decomposition (|13|) of the propagator in 
terms of the P, Q matrices which satisfy many interesting 
relations. We list below some of the relations that are 
useful for our discussion. 



QaaQ = 0, QasPiAn, Aa) = AaQ, P(Ar, AA)a3Q = ArQ, 

P(Ai,R, Ai,a)ct3P(A2,R, A2,a) = P(Ai,rA2,R, Ai,aA2,a). 



(19) 



r 



Using these relations, the n-point amplitude in H18|l can 
be simplified considerably. First, we note that the ex- 
pression on the right hand side can at most be linear in 
Q and, therefore, can only have at most a single propaga- 
tor of the type A-| which as we have seen can describe 

on-shell particles (see, for example, Furthermore, if 
we assume that an on-shell propagator can be thought of 



J 



as a cut open line (representing an on-shell particle), it 
is clear that the retarded n-point amplitude will involve 
only connected diagrams (not disconnected as can be the 
case in a time ordered Feynman amplitude which we have 
seen earlier). In fact, an explicit evaluation of (|18|l leads 
to 
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(20) 
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where we are identifying i„+i = ti,En+i = Ei. We are 
also using the convention that when m = (or m = 
n — 1), the term in the parenthesis has the value 



(21) 



It is obvious that the first term in the bracket that in- 
volves only a product of retarded propagators would van- 
ish when integrated. Furthermore, using (|12|l we can con- 
vert all the advanced propagators into retarded ones and 
write 
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(22) 
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This gives a forward scattering description for the re- 
tarded n-point amphtude at one loop at zero tempera- 
ture. Each term in the series is a number of retarded 

propagators with one on-shell propagator (A-| ) leading 

to the forward scattering of a single on-shell particle in 
all possible manner in a connected causal manner. Un- 
like the Feynman amplitude in the forward scatter- 
ing description for the retarded n-point function does not 
involve disconnected diagrams which is also reflected in 
the basic definition of the retarded amplitudes in terms 
of nested commutators that we will discuss in section III. 
From H22|l . we can easily derive a recursion relation for 
the integrands of the one loop retarded amplitudes of the 
form 

n 



i=l 



-h7i^^AR(-t„+i -Hii,£;„+i). 



(23) 



ward scattering description for some simple higher loop 
graphs can be easily derived algebraically as follows. Let 
us consider the scalar 0"+^ theory. In this case, the re- 
tarded self-energy at n-loops (see Fig. 5) can be written 
as 



\++{En+l, 



(2^ 
(2^ 
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.A+_(£;,)AR,(i;„+i). 



(24) 



B. Forward scattering amplitude at higher loops 

The simple algebraic derivation of the one loop forward 
scattering amplitudes for the zero temperature retarded 
n-point function does not carry over to higher loops in 
general. This is simply because of the fact that at higher 
loops, more than two propagators (internal lines) may be 
connected to a given vertex. In such a case, the conve- 
nient matrix structure for retarded amplitudes that arises 
in one loop (because only two propagators can be con- 
nected to a vertex) is not present. Nonetheless, the for- 



Here ki,i = 1,2, ■ ■ ■ ,n denote the n independent mo- 
menta of the loops and in the intermediate steps, we 
have used various relations such as and have ne- 
glected terms involving products of retarded quantities 
in the integrand (which will vanish upon integration). 
The recursion relation in H24() is interesting for two rea- 
sons. First, it shows the generic feature in higher loops 
that any retarded amplitude at n loops can be given a 
forward scattering description in terms of retarded am- 
plitudes at lower order. Second, the recursion relation 
can be thought of as a recipe for opening up loops 
|18|| for this particular diagram. 
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Although the forward scattering description for some 
simple higher loop diagrams can be derived algebraically, 
for a general higher loop amplitude, this is best estab- 
lished diagrammatically. For this purpose, let us intro- 
duce the graphical representation for the two parts of the 
matrix propagator in as 

P , — a b 

^ ah — 5 

A+_Qafc = ^ ^. (25) 

There are two important things to note here. First, the 
"cut" propagator corresponds to the on-shell propagator 
and all the elements of the matrix Q have the value unity 
so that the form of the "cut" propagator is the same in- 
dependent of the indices a, 6 = ±. Second, since the 
propagator P is directional, we choose the convention 
of taking the direction of time flow to be towards the 
"-I-" vertex in a retarded amplitude (which corresponds 
to the largest time). This simplifies the derivation and 
is physically meaningful to give a causal evolution for 
the amplitudes. The direction of the time flow at other 
vertices, where the thermal indices are being summed 



over, is unimportant. It is clear that a graphical decom- 
position of the propagator in the manner as described 
in H25|l would allow us to write any diagram as a sum 
of diagrams each consisting of certain numbers of "P" 
propagators and the remaining ones "cut" propagators. 
Each "cut" propagator corresponds to an on-shell prop- 
agator and, therefore, can be thought of as a cut open 
line representing forward scattering of an on-shell parti- 
cle. This can, therefore, also be thought of as a graphical 
description of the opening up of loops. 

In this process of "opening up of loops" , we may run 
into disconnected diagrams. As we will discuss in more 
detail in the next section, the retarded amplitudes corre- 
spond to vacuum expectation values of products of nested 
commutators and as such cannot have disconnected di- 
agrams (which would correspond to products of vacuum 
expectation values). This can, of course, be checked 
graph by graph at any order as was done at finite temper- 
ature in 'l9|. However, this can also be seen graphically 
as follows. Suppose, in this process of opening up of a di- 
agram, it separates into two disconnected parts as shown 
in Fig. 6. 



J 



<> o 

^ C 

C 

FIG. 6: Two ways of cutting the two-loop self-energy diagram with two cut propagators which render the graph disconnected. 
The thermal indices a, h and c are being summed 



In this case, there are two distinct possibilities. First, one 
of the disconnected parts contains the "-I-" vertex corre- 
sponding to the largest time and the other a connected 
part involving only vertices whose thermal indices are 
being summed over. In this case, the second part would 
vanish because of the identity l|16(l . The other possibil- 
ity is that one of the disconnected parts is a connected 
diagram involving the "-I-" vertex and the other simply 
consists of a vertex whose thermal index is being summed 
over. Once again, when we sum over the thermal index 
of this disconnected vertex (with a fixed distribution of 
the other thermal indices), the diagram would sum to 
zero (since the "4-" and the "— " vertices have a relative 
negative sign) . The crucial ingredient that allows this ar- 
gument to go through is the special property that a "cut" 
propagator is the same for any value of the thermal in- 
dices. As a consequence of this nice result, it follows now 



that for an arbitrary retarded amplitude at n loops, there 
can at the most be n number of "cut" propagators in a 
diagram because more cuts than that would render the 
graph disconnected. Furthermore, only those propaga- 
tors in a diagram can be "cut" propagators (even when 
their number is less than or equal to n) if they do not 
render the diagram disconnected. The propagators that 
are not cut correspond to the "P" propagators which can 
be seen explicitly from to be lower triangular with 

P++ = Ar, P+_ = 0, Pr = P++-P+- = Ar. (26) 

As a result, in a retarded amplitude, the uncut propaga- 
tors simply correspond to Ar propagators (which can be 
thought of as retarded "P" propagators). 

From these interesting properties follows a simple 
recipe for constructing the forward scattering amplitudes 
for a retarded graph at n loops. Start with a given graph 
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and write it as sum of all possible diagrams involving 
"uncut" and "cut" propagators (open lines) such that 
none of the diagrams is disconnected and that there are 
at the most n number of "cut" propagators. The diagram 
with n number of "cut" propagators would correspond to 
a tree level forward scattering diagram with intermedi- 
ate retarded propagators. Any diagram with the number 
of "cut" propagators less than n, would involve vertex 
diagrams of lower order (loop) as well as intermediate 
propagators that are retarded. The vertex diagrams of 
lower order would correspond to retarded diagrams with 



respect to the "P" propagators and, therefore, would in- 
volve only Ar propagators. This is the forward scatter- 
ing description for a retarded diagram at any loop at zero 
temperature. 

The above recipe is already obvious in the examples 
that we have discussed before. Let us illustrate these 
as well as some nontrivial examples at higher loops in 
a graphical manner. First, let us look at the one loop 
retarded self-energy in the <p^ theory which can be written 
as 



J 




+ R 



(27) 



Similarly, the retarded three point function at one loop 
in the (j)^ theory takes the form 



r(i) - 

^3,R ~ 




^ 3,R 



(28) 



Let us next look at the retarded self-energy at two loops 
in the (j)'^ theory, which takes the form 
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Here and in what follows a multiplicative factor in a at a nontrivial diagram for the retarded self-energy at two 
graph denotes symbolically the number of distinct graphs loops in the (j)^ theory which takes the form 
of the same topology that can be drawn. Let us next look 




_ y(2)(P) 




T 



{3).R 



R ^ R R 



(30) 



All these examples illustrtate how the recipe works for 
an arbitrary retarded amplitude at 71-loops and demon- 
strate the forward scattering description for a retarded 
amplitude at zero temperature. 



III. FORWARD SCATTERING DESCRIPTION 
AT FINITE TEMPERATURE 

Given the forward scattering description for retarded 
amplitudes at zero temperature, it is now straight for- 
ward to derive the forward scattering description at fi- 
nite temperature through the use of the thermal opera- 
tor. Let us recall that in the closed time path formalism, 
the thermal propagator for a scalar field can be related to 
the zero temperature one through the thermal operator 
as [ll 

A^^\t,E) = 0^^\E)A{t,E), (31) 

where 

O^^HE) = l + nBiE){l~SiE)). (32) 

Here nB{E) represents the bosonic distribution function 
and S{E) is a reflection operator that takes E —E. 
The thermal operator is the same for each component of 
the propagator (it is a scalar multiplicative operator) and 



leads to 

A^+l{t,E) = 0^^\E)A+^{t,E) 

= A+^{t,E) + Al_{t,E), (33) 

where we have identified the temperature dependent part 
of the propagator to be (this notation is in an attempt 
to be consistent with the notation in although we 
have denoted the propagator in those papers by G) 

A^_(i,£;) = ^(e-^* + e^^'). (34) 

The other interesting thing to note is that 

A^PiUE) = 0^{E)An.it,E)^Anit,E), 
A^Pit,E) 0^^\E)AAit,E) = AAit,E). (35) 

Namely, physical propagators such as the retarded and 
the advanced propagators (at the tree level) are indepen- 
dent of temperature. Graphically, these results can be 
written as 

0^^\E)^-^„ = 

0<^^HE),r^^„ = .r^-. + (36) 

where the "double cut" propagator can be identified with 
A^_Q (completely parallel with the notation of 0|). 
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We know that any Feynman graph at finite temper- 
ature is related to the corresponding zero temperature 
graph through a thermal operator 0, |2| that can be built 
out of the basic thermal operator in (|32|l . For example, 
the integrand of a graph (after the internal time inte- 
grations are done in the mixed space or the energy inte- 
grations are done in the energy-momentum space) with 
N scalar propagators carrying energy Ei,i = 1,2, N 
at finite temperature is related to the integrand of the 
corresponding graph at zero temperature by the thermal 
operator 

N 

O^^^ =Y{0(^\e{). (37) 

i=l 

This is a consequence of the simple fact that at finite tem- 
perature only the propagators of the theory are modified 
(because of the periodicity properties) while the inter- 
action vertices remain unaltered. As a result, the finite 
temperature forward scattering description for a retarded 
diagram can be obtained directly from the zero temper- 
ature one by simply applying the thermal operator ap- 
propriate to the particular diagram. Of course, this also 
clarifies the origin of the finite temperature forward scat- 
tering description, namely, it exists because there is a 
corresponding description at zero temperature. 

It is clear from that since the thermal operator 
does not change the "P" propagators, all the uncut lines 
in a diagram in the forward scattering description will 
continue to be the zero temperature retarded propagator, 
Ar. The thermal operator will only change the "cut" 
propagators (the open lines) to a "cut" plus a "double 
cut" propagator. (Like the "cut" propagators, the "dou- 
ble cut" propagators are also on-shell with a factor of 
the distribution function tib and correspondingly can be 
thought of as representing thermal on-shell incoming and 
outgoing particles.) Thus, one can organize the graphs in 
the number of "double cut" propagators. There will be 
diagrams with no "double cut" propagator, a single "dou- 
ble cut" propagator and so on and the maximum number 
of "double cut" propagators will be n for a retarded am- 
plitude at n loops. The diagrams without any "double 
cut" propagator will, of course, correspond to the zero 
temperature retarded amplitude. The diagrams with the 
maximum number of "double cut" propagators {n in the 
case under study) will represent tree level forward scat- 
tering amplitudes for thermal on-shell particles. The dia- 
grams with the number of "double cut" propagators less 
than n will all arrange into forward scattering amplitudes 
for thermal on-shell particles with (zero temperature) re- 



tarded vertices of lower order (n — 1 and lower) . Namely, 
the effect of applying the thermal operator to a retarded 
amplitude at zero temperature is to change all the "P" 
retarded vertices in the forward scattering description to 
genuine retarded vertices at zero temperature and replace 
all the on-shell forward scattering particles by thermal 
on-shell forward scattering particles. If we ignore the zero 
temperature retarded amplitude, the rest of the diagrams 
yield the temperature dependent forward scattering am- 
plitudes. 

The fact that the graphs will arrange as described 
above can be seen symbolically as follows. Let us identify 

Ar = P, A^ = y,A^_ = y^. Then, a retarded graph 

with TV propagators at zero temperature can be symbol- 
ically represented as (P + y)^ . If the graph is at n loops, 
then we can expand it in terms of the number of on-shell 
propagators {y) and write symbolically as 

^ = iP + yf=P'' + aiP''-'y+--- + a,,P^-y". 

(38) 

Here the multiplicities ai,i = 1,2, ■■■ ,n are assumed to 
denote the number of ways the forward scattering on- 
shell particles can occur in a graph without disconnecting 
the diagram (which is why these cannot correspond to 
the pure binomial coefficients). Furthermore, the terms 
involving powers of P represent intermediate retarded 
propagators as well "P" retarded vertices. As we have 
seen, under the action of the thermal operator 

0(^)(P + y) = (P + y + /). (39) 

As a result, using the thermal operator representation, 
we obtain 

p(»)(T) _ 0(T)p(n) - (p + y + yP]N 

= (P + y)^ + ai(P + 2;)^-i/ + ... 
+ --- + a„(P + y)~-"(/)" 

~ ^ (N),K^ y (Af-1),R + (7V-2),R + ) 

+(/)^WriH+---)+---- (40) 

It is useful to remember that the sum of the power of y 
and y^ in any term on the right hand side of H40|l can at 
the most be n at rt loops. 

Let us illustrate these relations with some examples. 
Applying the thermal operator to the forward scattering 
description of the retarded self-energy in the (j)^ theory 
at one loop (see (123), we obtain 
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(i)(P) 

R 



.(1) 



.(1) 



-(l)/5 



(41) 



I 

Similarly, the one loop retarded three point function (see l|28ll l at finite temperature takes the form 
I 



r. 



3,R 



R R , 



R R 



= r 



3,R 



R R , 



^ ^ 3,R ~ ^ 3,R 



(1) ^ \. r/ + \. + ^/ =r« 



(42) 



I 

A slightly more complicated example would correspond (I29II 1 at finite temperature which takes the form 
to the two loop retarded self-energy in the 0^ theory (see 



.(2)(T) _ 
^R — 




(2)/3 



(43) 



This shows explicitly how the "P" retarded vertices of 
lower order rearrange themselves into full retarded ver- 
tices. Finally, let us consider the nontrivial example of 



r 



the two loop retarded self-energy diagram (see 1(201)) 
finite temperature 
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This nontrivial example once again demonstrates how the 
"P" retarded vertices of lower order rearrange themselves 
into full retarded vertices. We note that the thermal op- 
erator C'(^) for the different amplitudes in iflTl . . 
and (|44|l are different and their appropriate forms can be 
obtained from (|37|l . Furthermore, the temperature de- 
pendent forward scattering amplitudes are denoted with 
a superscript (3 to coincide with the definition in |l9j | and 
in the applicable examples can be checked to agree with 
the results there. 

We would now like to make some observations on the 
structure of retarded amplitudes in general which are not 
directly related to the main goal of this paper, but are 
quite important in understanding their structures. First, 
we note that the retarded iV-point amplitude at any loop 
is defined algebraically in the coordinate space in terms 
of the original fields of the theory as the vacuum expec- 
tation value of the nested commutators |23 

rAr,R(tl, ^2, • • • , tN) 
= {~lf-^0{h - t2)9{h - is) • • • 0{tN-l - tN) 

x(O|[[---[</.(a:i),(/.(x2)],0(a:3)]--- ,0(xAr)]|O) 
-|-permutations. (45) 

Here we have assumed that the time coordinate ti is 
the largest among all the coordinates and the "permuta- 
tions" refer to symmetrizing in all the other coordinates 
(other than the largest time) and fields. As a result, the 
retarded amplitude is symmetric in all the coordinates 
other than the largest time coordinate. For a real scalar 
field, (t>{x) is a Hermitian operator and, therefore, the 



factor (— i)^^^ shows that the retarded amplitudes are 
real in coordinate space. Under Hermitian conjugation, 
the change in the sign of each factor of i is compensated 
by the change in sign coming from each commutator. We 
have already argued in section II graphically that a re- 
tarded amplitude cannot have disconnected parts which 
would correspond to products of vacuum expectation val- 
ues. This can also be seen from the above definition as 
follows. Let us denote the nested commutator involving 
the first — 1 fields as 

A - [[• • • [0(xi), 0(X2)], (/.(Xg)] • • ■ , ^{XN^I)] . (46) 

Then, we can write the vacuum expectation value of the 
nested commutator in (|45|l as 

(0|[A,(/)(xAr)]|0). (47) 

Inserting a complete set of intermediate states (say, dis- 
crete energy eigenstates or the particle number states), 
this can be written as 

{0\[A,^{xn)]\0) 

OG 

= ^ mA\n){n\^{xN)\0) - m{xN)\n){n\A\0)) 

OG 

= J2 {m\n){n\^{xN)\0) - m{xN)\n){n\A\0)) . 

(48) 

Namely, the intermediate vacuum states cancel out in 
the vacuum expectation value of the commutator and as 
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a result, the retarded amplitudes contain only connected 
graphs which we have explicitly seen earlier. 

It was noted earlier that the retarded self-energy for 
a real scalar field in the mixed space is a real quantity. 
As we have already argued, the retarded amplitudes, by 
definition, are real in the coordinate space. However, in 
going to the mixed space, one Fourier transforms the spa- 
tial coordinates into spatial momenta and Fourier trans- 
formation does not maintain the reality of a function in 
general. Let us comment here briefly on when the re- 
tarded amplitudes will be real in the mixed space for a 
scalar theory. The definition of the retarded amplitudes 
in the mixed space take the form (see MSIl ) 

^N,R{tl,t2, • • • , tj^i) 

= {-if-^e{ti - h)e{t2 - is) • • • 0{tN-i - tN) 
x(o| [■ • • m^pi), 0(^2, P2)], • • • , HtN^m)] |0) 

-l-permutations. (49) 

Under Hermitian conjugation, the change in sign in each 
factor of "i" is still compensated for by the change in sign 
coming from each commutator. However, since under 
Hermitian conjugation 

0(i,p)->0(i,-p), (50) 

the retarded amplitude is not real in general. In fact, let 
us note explicitly that under Hermitian conjugation, 

= i-i)^-'9{t, - t2)9it2 - is) • • • 0{tN-l - t^) 

X(0| [■ • • [0(<i, -pi), (/)(t2, -P2)], • • ■ , HtN, -Pn)] |0) 

+permutations. (51) 

We note that if the scalar field transforms under parity 
as 

^{t,p) -> 0(t, -p) = (-l)«0(t,p-), (52) 

where 77 = (—1)" denotes the intrinsic parity of the scalar 
field, then we can write 151(1 as 

rAr,R(*ij *2, • • • , tw) 

= (-1)^"(-Z)^-I0(tl - t2)e{t2 - is) • • • 0{tN-l - tN) 
X(0| [■ • • Pi), q>it2,P2)], • • • , 0(ijv,Piv)] |0) 

-|-permutations 

= (-l)^"fAr^R(il,i2,--- ,ijv). (53) 

For a scalar field of even parity, a — and we see that the 
retarded amplitudes will continue to be real even in the 
mixed space. However, for a pseudoscalar field, a = 1 
and we note that only parity conserving retarded am- 
plitudes will be real while the parity violating retarded 
amplitudes will be purely imaginary in the mixed space. 
We would like to emphasize that the reality of an am- 
plitude in the coordinate space/mixed space is not con- 
tradictory to the existence of dispersion relations in the 



energy-momentum space since the imaginary parts of the 
amplitudes in energy-momentum space arise from the 
imaginary parts of the step functions {0{t)) in the in- 
tegral representation. 



IV. FORWARD SCATTERING DESCRIPTION 
FOR YANG-MILLS THEORY 

The results of the earlier sections show that the for- 
ward scattering description at finite temperature can be 
obtained from the forward scattering description at zero 
temperature by the use of the thermal operator. Al- 
though we have done this explicitly for scalar field the- 
ories, this can be generalized easily to other theories. 
We would like to emphasize that this correspondence be- 
tween the finite temperature and zero temperature for- 
ward scattering descriptions should not be thought of 
as useful only in establishing a graphical identification. 
It is also quite useful as a calculational tool as well as in 
clarifying various aspects of field theories at high temper- 
ature. To give an example of this, we will next derive the 
retarded gluon self-energy in the Yang-Mills theory (be- 
longing to SU{N)) at one loop in the hard thermal loop 
approximation from the zero temperature result, which 
will also clarify the structure of this thermal amplitude. 

A lot is known about the structure of Yang-Mills the- 
ories at high temperature JT„ J^. It is known, for ex- 
ample, that in the hard thermal loop approximation the 
one loop retarded self-energy in the forward scattering 
description is independent of the gauge fixing parameter 
and has a manifestly gauge invariant (transverse) and 
Lorentz covariant structure (before carrying out the an- 
gular integrations). However, the reason for such a struc- 
ture at finite temperature is not well understood. We will 
see below that such a structure of the integrand for the 
retarded self-energy already exists at zero temperature 
in the appropriate regime and since the thermal opera- 
tor is gauge invariant, the thermal amplitude obtained 
through the application of the thermal operator repre- 
sentation preserves these properties. 

From the discussions in section II, we can immedi- 
ately write down the forward scattering description for 
the one loop retarded self-energy for the Yang-Mills field 
(including the ghost contributions) easily. Our discus- 
sion in section II has been completely within the context 
of the mixed space simply because we wanted to bring 
out the physical nature of the retarded amplitudes as 
evolving forward in time in a connected manner. How- 
ever, the thermal operator representation holds equally 
well in the energy-momentum space (which can be seen 
simply by Fourier transforming the external time coordi- 
nates) where the thermal operator acts on the integrand 
after all the energy integrations have been carried out. 
Since all the calculations of thermal amplitudes in the 
hard thermal loop approximation have been carried out 
in the energy-momentum space, in this section we will 
also work in the energy-momentum space and use the 
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thermal operator representation in this space. 

The forward scattering description for the one loop re- 



tarded self-energy for the Yang- Mills field (belonging to 
SU{N)) at zero temperature can be written as 



J 




P ('X R !(' ,o 
k'. c k+p, d k.c 
52 f d^k 



{27r)3 



k.c k+p.d k,c 



d^^{k+p) 

[k+pY 



-k) 



dP''{k)V,„^{~p,~k,k+p) 



I 



(54) 



where wavy and dotted lines in the graphs denote respec- 
tively the gluon and the ghost propagators, which in an 
arbitrary covariant gauge have the forms 



D''\k) 



df_t^ik) 



(55) 



Here ^ denotes the gauge fixing parameter and the appro- 
priate "ie" factors in the denominators of the propagators 
are understood. We also note here that a product such 
as 5+(fc2)p- in ifSUl has to be understood in a regularized 
sense as has been described in Similarly, the gauge 
and the ghost vertices have the forms (we are suppressing 
the energy-momentum conserving delta functions and are 
gr'^V^uX and V^"= = gr'^V^) 



defining V^^^ = 



Vp^^\{p,k,q) = i]f^iyip-k)x + rii^\ik-q)^ + rix^{q-p)^, 



Vf,{p,k,q) = fc^. 

The quartic vertex, {-ig^W^lf^) 
l22ll. Furthermore, we have defined 



(56) 

can be read out from 



(5+(fc2) =6i(fco)(5(fc^). 



(57) 



Since the ghost propagator as well as the gauge interac- 
tion vertex are independent of ^, for the purpose of un- 
derstanding the ^ dependence in H54|) we can restrict our- 
selves only to the diagrams involving intermediate gauge 
field propagators. 

The tree level gauge field vertices satisfy the iden- 
tity 

k''V^,y\(J), k, q) = (jy^AP^ - Pf,P\ - t^^a'?^ + ^m^a) , (58) 

where we have used the fact that k = —q — p because 
of energy-momentum conservation. In the region where 
kfi 3> p^ (where we are assuming that k, q are the internal 



momenta and p is the external one) , the first two terms in 
(|58|l can be neglected for our purpose and what remains 
is a structure that is transverse to the other internal mo- 
mentum. As a result, we see that the terms quadratic 
in the parameter ^ vanish so that the self-energy in H54|) 
can at most depend linearly on the gauge fixing param- 
eter ^. An explicit evaluation of the diagrams involving 
the gauge field propagators shows that the linear terms 
in ^ cancel out among the two classes of diagrams at the 
integrand level. Consequently, in the limit fc^ ^ P^j., the 
leading order contribution to the self-energy is indepen- 
dent of the gauge fixing parameter ^. 

At zero temperature, the amplitudes are manifestly 
Lorentz covariant. Furthermore, the Ward identity for 
the gluon self-energy requires that it be transverse to the 
external momentum, 



p'K}^2iP) = 0- 



(59) 



Since the leading part of the retarded self-energy is in- 
dependent of the gauge fixing parameter ^, in evaluating 
(|54|l in this regime, we can choose it to have the value 
^ = 1 for simplicity (Feynman gauge). In this case, the 
numerator of the gauge propagator in H55|) is simply the 
metric tensor independent of any momentum. The mo- 
mentum dependence in the numerator in (|54|l comes from 
the vertices and is, therefore, at most quadratic in the 
momenta. (Note that the three point vertex depends on 
the momentum, but the quartic vertex is independent of 
the momentum.) Since we are interested in the "hard" 
internal momentum region where k^ ^ p^, we can ex- 
pand the numerator in powers of the internal momentum. 
Similarly, we can also expand the denominator in (|54l) in 
this region as (recall that k-^ — because of the delta 
function) 



(60) 



{p + fc)2 2p ■ k {2k- pY 



Using these in H54|l . one explicitly finds that the lead- 
ing contribution in this region comes from terms in the 
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integrand (multiplying the delta function) which are of 
degree zero in both the external as well as the inter- 
nal momentum. There is only one such tensor structure 
available at zero temperature which is consistent with the 
Ward identity (|5^ . namely, 



P 



k 



p^k^ ki, 

(p ■ fc)2 ■ 



(61) 



and the integrand has to be proportional to this struc- 
ture. Explicit calculation indeed shows that in this re- 
gion, H54|) takes the form 



Pfj.ki/ - 



Pi^kfj^ p kp^ky 



(62) 



This is the leading contribution in the "hard" internal 
momentum region and as we will show shortly, this leads 
to the hard thermal loop results at high temperature. 

There are several things to note from the structure 
in H62|l . First, the quantity in the parenthesis is mani- 
festly Lorentz covariant of degree zero in both the inter- 
nal and the external momenta. Furthermore, it is man- 
ifestly transverse (gauge invariant). We note that the 
integral in (|62|l is quadratically divergent and is usually 
set to zero in the dimensional regularization. However, 
if we only carry out the fep integration it does not van- 



ish and as we will show shortly, it is this leading con- 
tribution that leads to the hard thermal loop results at 
high temperature. (We remark here parenthetically that 
the conventional ultraviolet divergent terms in the self- 
energy are, in contrast, gauge dependent and are only 
logarithmically divergent, yiel ding at high temperature 
to only InT contributions y^.) We want to emphasize 
that the reason for carrying out only the fco integration 
is that the thermal operator acts on the integrand after 
the energy integrations have been carried out and before 
evaluating the integrations over the spatial momenta . 
To apply the thermal operator representation, we need 
to integrate over the fcg variable in l|62() and this leads to 



2Ek 



(27r)3 

Pfj,kiy -\- Pvk^ P^k^ki/ 



p ■ k 



{p-kf 



(63) 



where — \k\ and we have defined fc^ — {1,—k). Apply- 
ing now the thermal operator to the integrand, we obtain, 
(It is worth clarifying here once again that although some 
of the diagrams in 1)54(1 involve two propagators, only the 
thermal operator corresponding to the on-shell propaga- 
tor leads to a nontrivial result. The retarded propagators 
are unchanged by the application of the thermal operator 
as discussed in (|36ll .') 



J 



n' 



ab{l){T) 



ah 



_ i-pab(l) 



ah 



d^k nsjEk) 
(27r)3 Ek 



Pt^k 



V 



p ■ k 



P^k^ ki, 



p ■ k 



[P ■ kf J 

Puk^ ^ p^k^ki, 
{p-k)\ 



n 



a6(l)/3 



(64) 



The temperature independent part n°^^p can now be set 
to zero using dimensional regularization for the spatial 
momentum. 

We note that since Ep; = \k\, the radial momentum 
integration in the second term in (|64|1 can be done to 
yield 



/>oo 

/ dk kn^{k) = 
^0 



6 



(65) 



which leads us to the temperature dependent part of the 
self-energy 



n 



dn 



487r 

P^ki/ -\- Pvk^ P^k^ky 



p ■ k 



{p ■ kf 



(66) 



This is the well known result for the temperature depen- 
dent part of the retarded self-energy in the hard thermal 
loop approximation in the forward scattering description. 
The angular integrations break the manifest Lorentz in- 
variance. (At finite temperature, Lorentz invariance is 
broken because the rest frame of the heat bath defines a 
preferred reference frame.) However, the beautiful prop- 
erties of the integrand follow simply from the properties 
of the zero temperature amplitude. 

This example demonstrates that, in addition to estab- 
lishing a graphical correspondence, the thermal operator 
representation can also be conveniently used to calculate 
the thermal amplitudes in the hard thermal loop approx- 
imation from the "hard" internal momentum amplitudes 
at zero temperature. Such a calculation also clarifies var- 
ious important features of the thermal amplitudes at high 
temperature. 
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V. CONCLUSION 

In this paper, we have derived systematicaUy the for- 
ward scattering description for retarded ampHtudes to 
all orders at zero temperature. This graphical derivation 
then allows us to obtain the forward scattering descrip- 
tion for such amplitudes to all orders at finite temper- 
ature through the thermal operator representation. Al- 
though our derivation has been within the context of a 
scalar field theory, the derivation can be generalized eas- 
ily to other theories with or without a chemical potential. 
Furthermore, although we have used the real time formal- 
ism of closed time path for our discussions for simplicity, 
the results also hold in the imaginary time formalism 
(which we do not go into). Besides giving a graphical 
derivation of the forward scattering description at finite 
temperature, such a relation can be used as a power- 
ful tool for calculations at high temperature and clarifies 
various properties of thermal amplitudes. As an exam- 



ple, we have calculated the one loop retarded self-energy 
for gluons in the Yang-Mills theory at finite tempera- 
ture starting from the forward scattering description at 
zero temperature. This derivation emphasizes that var- 
ious nice features of these amplitudes such as gauge in- 
variance, transversality, manifest Lorcntz covariance etc 
arise simply because the zero temperature amplitude al- 
ready possesses such properties. This description of the 
forward scattering amplitudes at finite temperature pro- 
vides yet another example of the usefulness of the thermal 
operator representation. 
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